Noise correlations of a strongly attractive spin-1/2 Fermi gas in an optical lattice 
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We calculate density-density correlations of an expanding gas of strongly attractive ultra-cold spin- 
1/2 fermions in an optical lattice. The phase diagram of the tightly bound fermion pairs exhibits 
a Bose-Einstein condensed state and a Mott insulating state with a single molecule per lattice site. 
We study the effects of quantum fluctuations on the correlations in both phases and show that they 
are especially important in the Bose-Einstein condensate state, leading to the appearance of singular 
peaks. In the Mott insulating state the correlations are characterized by sharp dips. This can be 
utilized in experiments to distinguish between these two phases. 



I. INTRODUCTION 

Experiments with ultra-cold atoms have opened many 
directions to study ground state properties of complex 
many-body systems, such as Bose-Einstein condensation 
and fermion pairing. The introduction of optical lattices 
has added the possibility of creating new ground states 
[1]. A prominent example is the formation of a Mott 
insulating (MI) phase for strongly interacting bosons [2] 
or fermions [3]. 

Another interesting direction was brought to the field 
by the unprecedented control over the interaction be- 
tween atoms via a magnetic Feshbach resonance [4]. By 
means of the latter it became possible to turn fermions 
with two hyperfine internal states into bosons [5] by pair- 
ing. For weak attraction the Bardeen-Cooper-Schricffcr 
(BCS) theory predicts an instability towards formation 
of large Cooper pairs, whose size can be much larger 
than the separation between atoms in the Fermi gas. If 
the attraction is turned on stronger, then these Cooper 
pairs become smaller and eventually, for very strong in- 
teraction, they behave as point-like (hard-core) bosonic 
molecules. Although there is no sharp distinction be- 
tween the BCS state and a condensate of local pairs of 
fermions, the cross-over from the weakly interacting to 
the strongly interacting Fermi gas is of substantial inter- 
est [6]. The experimental creation of bosonic molecules 
in a Fermi gas atoms as well as their Bose-Einstein con- 
densation has been recently reported [7] . This BCS-BEC 
crossover has also been realized in an optical lattice [8]. 
In the presence of an optical lattice the tightly bound 
molecules in the BEC regime may undergo a quantum 
phase transition from a BEC state to a MI state with 
one particle per lattice site. This leaves us with question 
how to identify the two phases by simple measurements. 

A clear signature of superfluidity in a system of bosons 
is the presence of sharp density peaks due to phase coher- 
ence in time-of- flight experiments [2], In the MI phase, 
on the other hand, the density peaks are washed out due 
to the absence of phase coherence (see [9]). To character- 
ize this phase it is useful to study the noise correlations 
in the Fermi gas [10], where a characteristic feature is a 
dip that has been seen experimentally [11]. 

In this paper, we study density-density correlations of 
an expanding gas of ultra-cold spin-1/2 fermions realized 



in an optical lattice. We discuss that the density-density 
correlations reveal the fermionic nature of the gas in the 
MI state and its bosonic nature in the Bose-Einstein con- 
densate (BEC). We show that the effects of quantum 
fluctuations on the correlations are important. 



II. MODEL 

An effective Hamiltonian for fermions in an optical 
lattice near a Feshbach resonance was derived in Ref. 
[12]. Under typical experimental conditions, the result- 
ing Hamiltonian can be mapped to an effective single- 
band model, which comprises the local dressed bosonic 
molecules and individual fermionic atoms separately. In 
this paper we consider a purely fermion model in an op- 
tical lattice with attractive interaction and with two in- 
ternal states, denoted by the pseudospin projections f 
and |, respectively. The fermionic field operator Vv(r) is 
expanded in a basis of Wannicr functions [13] as 



Vv(r) = ^ w„(r - Ri)ci 



(1) 



which allows us to use the tight-binding representation. 
Here, Ci >n ^ is the annihilation operator for particles in 
the Wannier state n at lattice site R^. For a deep optical 
lattice only the lowest Wannier state is taken into account 
with ci M = Ci.o.tr- In this case the lattice Hamiltonian 
reads 
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Nearest-neighbor tunneling of the individual fermions is 
described by the parameter i. Moreover, there is also 
a nearest-neighbor attraction with parameter J which 
can be understood as a tunneling term of fermionic pairs 
(molecular tunneling). A similar model of a mixture of 
hard-core bosons and fermions in a lattice was discussed 
in Ref. [14]. However, in our model the hard-core bosons 
hi are given by pairs of fermion operators CjfCij_. This 
means that the fermions can transmute dynamically into 
hard-core bosons by local pairing. 
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III. PHASE DIAGRAM 

We apply the functional integral formalism developed 
for the Hamiltonian in Eq. (2) in [15, 16]. The partition 
function of a grand-canonical ensemble of fermions with 
chemical potential [i at inverse temperature (3 can be ex- 
pressed in terms of an integral over a fermionic (Grass- 
mann) field ip and complex molecular fields <fi, \ as 

Z = J e-SS^^x] (3) 

with the Lagrangian 

c = E + ^7 E xiXi E > (4) 

id i id 

where we have used the notation & = (ijj^, ipi) and C = 
(ipl, "0|) T , and with the inverse fermionic Green's matrix 

6 _ 1= / -Ufi-X -^It*). ( 5 ) 

\d T -fi-t i<j) + x J 

Here Vij = 4g&\i— t is the nearest-neighbor tunneling 
matrix with elements i/2d. 

The integration over the complex molecular fields (j>, x 
can be performed in saddle-point approximation [15, 16]. 
For a vanishing fermionic tunneling, i.e. by assuming 
a situation where all fermions are paired up, the phase 
diagram is depicted in Fig. 1. There are three phases: 
the BEC of molecules with the condensed fraction uq — 
(J 2 — 4/i 2 )/(4J 2 ), the MI state with one particle per site 
and no = and the empty phase. Fluctuations around 
the saddle point provide the low-energy excitations of the 
bosonic molecules. They are gapless in the BEC phase: 

£k = \J^J 2 9^n +Aglu 2 , (6) 

where <?k = 1 — l/d^i—i cosfc^ is the dispersion of the 
free Bose gas. On the other hand, the excitations of the 
MI state have a gap A = 2u — J > : 

e k = A + J 9k . (7) 



IV. DENSITY-DENSITY CORRELATIONS 

In a typical experiment with ultra-cold atoms it is not 
easy to identify the spin projection. Therefore, the de- 
tected fermion density should be described by a superpo- 
sition of both spin projections. For this purpose we intro- 
duce the spin-independent fermionic density n = n-f+rij. 
Then the truncated density-density correlation function 
in a time-of-fiight experiment is given under the assump- 
tion that the atoms expand freely after they have been 
released from the optical lattice as [17] 

= (n(r)n(r'))-(n(r))(n(r')> - 




FIG. 1: Phase diagram for t = and k B T = 0. The three 
phases are the BEC with a non- vanishing condensed density, 
the MI states with vanishing condensed density and with one 
molecules per lattice site and the empty phase, n and J are 
in arbitrary energy units. 

f-J (n(k)n(k'))-(n(r))(n(r')). (8) 

Here k is related to r by k = Mr j fit with the effective 
mass M/h oc 1/J. Free expansion assumes that atoms 
evolve independently. This can be achieved by switching 
off an optical lattice as well as by switching the magnetic 
field to values far from the strongly attractive regime 
suddenly before the expansion [10, 18]. The above for- 
mula reveals then that the density distribution, observed 
after a fixed time-of- flight, is related to the momentum 
distribution of the particles trapped in the lattice before 
the expansion. It is widely used to analyze ground state 
properties of ultra-cold atoms trapped by an optical lat- 
tice [19]. 

Now we use the expansion of Eq. (1) to express the 
density-density correlation of Eq. (8) in terms of the 
coordinates of the underlying optical lattice as 

Cr '" = (|Q 2 V( k )l 2 H k ')l 2 x 

i,i',jd' 

+(n(r))5(r-r')-(n(r))(n(r')) (9) 

where w is a Fourier transform of the Wannier func- 
tion. Using the saddle-point approximation of the com- 
plex fields (4>,x) we can evaluate {c[ a c'jgCj'pCii a ). Since 
the fermion field appears only in a quadratic form in Eq. 
(4), the integration are given by a Wick contraction of 
pairs of fermions: 

{c\ a c) p C r i3C tla ) = {c] p Cj>l3) ~ (cl a C r0 )(cl (3 C l ' a ) + 
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(10) 



The fermionic expectation values on the right-hand side 
are matrix elements of the fermionic Green's function G 
ofEq. (5): 



and 



(12) 

where we have summed over Matsubara frequencies tu n . 
These expressions include also quantum fluctuations of 
the bosonic molecules when we take into account Gaus- 
sian fluctuations around the saddle-point solutions. For 
instance, we can integrate over these fluctuations to eval- 
uate the correlation function 



In general, the Green's function Gij is not diagonal 
in space due to the presence of the fermionic tunneling 
i but the decay is exponential with distance |R, — TLj\. 
The exponential behavior can be approximated by taking 
only the diagonal part into account. In other words, in a 
strongly interacting Fermi gas the formation of bosonic 
molecules and their dynamics is the dominant feature, 
and the individual tunneling of fermions is negligible. 
Moreover, the Green's function depends on the fields 
and x- I n that case the density-density correlations re- 
duces to 



2,/ 



Rk)| 2 R0)| 2 x 



RESULTS 



The expressions (cj|Cj|) and (c^Cii) arc constant due 



to translational invariance: 



C r o oc N 2 S. 



-1 



N AT 
N 2 



■5(k), (16) 



where we have denoted 2N = 2^-(c^Cj|) as the to- 
tal number of fermions and NgAf = 53i(^»T^*i) I 2 oc 
Af 2 \i(f> + x\ 2 ■ -A/" is the number of lattice sites. The ratio 
no = Nq/N is the condensate fraction, i.e. the relative 
contribution of condensed molecules. 

BEC.- For momenta close to the reciprocal lattice vec- 
tors G the main contribution in the BEC comes from the 
term <S(k), since it is singular for ~ 0: 



C r 



o cx 



5(k) 



4 Jn + Jgk 



(1 - 2 5k ) . (17) 



This results connects the density-density correlation 
function with the spectrum of the molecular condensate 
given in Eq. (6). The infrared divergence of the momen- 
tum distribution is a general property of a BEC at T = 
[20]. The related 2D plot of the correlation function is 
shown in Fig. 2. 

MI.- For larger densities the phase coherence in the 
molecular state is destroyed and a MI state with one 
bosonic molecule (i.e. a pair of fermions) appears. Due 
to strong phase fluctuations the second term in Eq. (16) 
vanishes, and the correlation function becomes 



C r ,0 OC 



k.G 



■5(k) 



(18) 



,ikR, 



E 



e* kR *(c iTCii ) 



S(k) 



(14) 

The first two terms constitute the density-density correla- 
tion function for fermions [17], whereas the third and the 
fourth terms account for the presence of the condensed 
molecules. In particular, the fourth term describes the 
effect of quantum fluctuations and can be expressed as a 
momentum distribution of the molecules [16] 



S(k) = £' 



k(Ri-Rj 



<W;» (15) 



with b\ = ct,cj|. Phase fluctuations can destroy these 
terms, for instance, in the case of a MI state of the 
molecules. It is important to notice that fermionic terms 
contribute with a negative sign in Eq. (14), in con- 
trast to the phase-coherent molecules, which contribute 
with a positive sign. This indicates a competition of the 
fermionic and the molecular contribution to the density- 
density correlation function. This provides a concept for 
measuring the properties of a strongly interacting Fermi 
gas. 



with the non-singular term 



J 2 
4/x 2 



5(k)«^(l-2 flk ). 



(19) 



Since the number of fermions N is large in Eq. (18), 
the contribution of S(k) is negligible in time-of- flight ex- 
periments of the MI. Thus the MI state is characterized 
by sharp dips in contrast to the singular peaks in the 
BEC, which appear on positions corresponding to the 
reciprocal lattice vectors G, similar to a Bragg diffrac- 
tion pattern. However, in contrast to the latter, where 
the diffraction pattern is created by light scattering on 
atoms of a crystal, here the pattern is created by the 
atomic state itself in the expansion process of the cloud. 

Eq. (19) describes the bosonic nature of the molecules 
in the MI state: There are hole excitations in the MI 
phase (particle excitations are suppressed by the hard- 
core nature of the bosons) and their dynamics contributes 
to the interference pattern via the factor oc f/k- A simi- 
lar situation was observed experimentally in the case of 
bosonic Mott insulator in Rcf. [21]. The 2D plot of the 
correlation function is shown in Fig. 3. 
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FIG. 2: Density-density correlation function of a strongly at- 
tractive spin- 1/2 Fermi gas in an optical lattice with J = 1, 
(i — 0.4, fcsT = 0, where the gas forms a BEC. The singu- 
lar behavior at the positions corresponding to the reciprocal 
lattice vectors is due to long range phase coherence of the 
condensed molecules. There is also small modulation due to 
the tunneling of molecules. The axes are given in units of 
ht/Ma, where a is a lattice spacing. 
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FIG. 3: Density-density correlation function as in Fig. 2 but 
with J = 1, /x = 0.6, ksT = 0. The singular behavior in Fig. 
2 is replaced by dips here. This is due to the anticommuting 
nature of fermionic operators and thus reflects the fermionic 
nature of the MI state. 




VI. DISCUSSION 



The density-density correlation function of an expand- 
ing cloud provides a clear picture of the state when it 
was still trapped by an optical lattice. It consists of four 
different contributions in Eq. (14), two of them are re- 
lated to the fermionic nature of the atoms. The second 
term leads to a fermionic dip in the density-density cor- 
relation function in the MI state. A third term mea- 
sures the condensate fraction in the case of condensed 
bosonic molecules and together with the singular fourth 
term leads to the sharp peaks. Even though the third 
term can compete with the second term, the singular 
term is the most relevant in the BEC phase. 



The competing behavior of the fermionic dips and the 
condensate peaks in the density-density correlation func- 
tion are a result of the anticommuting properties of the 
fermionic operators. This provides a simple concept to 
distinguish different states in a cloud of attractively in- 
teracting fermions. This behavior is rather different in a 
bosonic cloud, where all atoms contribute with the same 
sign to density-density correlation function because the 
bosonic operators commute [22]. 



In conclusion, we have studied an expanding cloud 
of strongly interacting spin-1/2 fermions after its re- 
lease from an optical lattice. The properties are de- 
scribed in terms of the density-density correlation func- 
tion. At lower densities a BEC is formed by the paired 
fermions, visible in the density-density correlation func- 
tion as sharp peaks. At higher densities a Mott-insulating 
phase appears, characterized by dips in the density- 
density correlation function. This distinct behavior can 
be used in future experiments to distinguish between 
these two phases in a gas of strongly attractive spin-1/2 
fermions in an optical lattice. 
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